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We demonstrate that a new tool, a model independent numerical Eliashberg inversion of the op- 
tical self-energy, based on maximum entropy considerations can be used to extract the magnetic 
excitation spectra of high-transition-temperature superconductors. In Bi2Sr2CaCu2 08+a we ex- 
plicitly show that the magnetic mode that dominates the self-energy at low temperatures directly 
evolves out of a smooth transfer of spectral weight to the mode from the continuum just above 
it. This redistribution starts already at 200 K in optimally doped materials but is much weaker in 
overdoped samples. This provides evidence for the magnetic origin of the superconductivity and 
presents a challenge to theories of the spin susceptibility and to neutron scattering experiments in 
high-temperature superconductors. 

PACS numbers: 74.25.Gz, 74.62. Dh, 74.72.Hs 



The phenomenon of high-temperature superconductiv- 
ity, discovered in the copper oxides twenty years ago , 
continues to challenge both theorists and experimental- 
ists. From the beginning it was clear that magnetism 
was an important part of the solution to this puzzle. 
In conventional superconductors the spectrum of pair- 
ing excitation was successfully extracted as an electron- 
phonon spectral function using an inversion of experi- 
mental tunnelling Q and optical data based on the 
Eliashberg equation. In principle this spectrum con- 
tains the information on the microscopic interactions 
among electrons mediated by boson exchanges which are 
needed to describe the superconductivity and can be cal- 
culated from band structure information Q. The pri- 
mary tools used to map out the magnetic excitations 
in high-temperature superconductors have been neutron 
scattering and nuclear magnetic resonance. The pic- 
ture that has emerged from these studies is a not-well- 
understood spectrum of excitation dominated by a con- 
tinuous background extending to unusually high ener- 
gies, sometimes called the Millis-Monien-Pines (MMP) 
spectrum H, IS 0> IE 01 which evolves into a broad peak 
in the local (q integrated) susceptibility at low tempera- 
ture along with a sharp resonance in the superconducting 
state centered at q = [w, tt), the 41 meV mode, named 
after its frequency in optimally doped YBa2Cu306+a 
(Y-123) [3 13 13 13 13 13 This sharp reso- 
nance has also been seen in T E Ba ^ CuOfi^A and in 
BiaSraCaCu^CW (Bi-2212) Hi 03. The magnetic res- 
onance [S 13 13 EJ E2 mode-phonons [H El Ef| 
and the broad continuum 0, 13 13 E3 m the bosonic 
spectra have been proposed as candidates for the pair- 
ing glue in the copper oxides. Here we use a broader 
definition of the magnetic resonance as the peak that de- 
velops in the local magnetic susceptibility at low temper- 
ature. Some theories of the spin susceptibility have pre- 
dicted that the magnetic resonance develops out of the 



MMP continuum [S |3| when the superconducting gap 
forms. However, so far no consensus has been achieved. 
To throw new light on this issue we have undertaken a 
study of the spectrum of excitations responsible for the 
self-energy of the carriers as a function of temperature 
and doping in three Bi-2212 systems using optical spec- 
troscopy focussing on the overdoped region of the phase 
diagram. We clearly observe experimentally for the first 
time the evolution of a strong peak in the integrated mag- 
netic response developing from the MMP continuum just 
above it as a function of both temperature and doping. 

The optical self-energy, Yj° p {lo) which involves a mo- 
mentum average, is closely related to the quasiparticle 
self-energy measured by angle resolved photoemission 
spectroscopy which is momentum specific. The optical 
self-ener gy is is defined in terms of an extended Drudc 
formula 1261: 



o-(w) 



(1) 



where io v is the plasma frequency and YP v (uS) = E^ p (a;) + 
iE 2 (u>). The optical scattering rate is defined by 
1/t° p {lo) = — 2Y,2 P (u)- The optical conductivity, a(ui) 
can be found through Kramers-Kronig transformation of 
the reflectance which is measured directly. The numeri- 
cal inversion of the optical scattering rate is based on an 
Eliashberg formalism. We start with a deconvolution of 
the approximate relation |3(J, |3lJ , 



t°p(uj;T) 



dQK(u,n;T)I 2 X (^), 



(2) 



where T is temperature, K(u),Cl;T) is a kernel deter- 
mined from theory, and / 2 x(f2) is the bosonic spectrum. 
For the deconvolution we utilized the maximum entropy 
method, originated by E. T. Jaynes [3]. For the kernel 
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FIG. 1: Optical scattering rate, l/r(w) as a function of u). 
Solid lines are data and dashed lines the Eliashberg fit. Heavy 
line is for temperature of 300 K (normal) and light for 27 K 
(superconducting) for optimally doped Bi-2212. In the inset 
we display the optical scattering rate of our two overdoped 
samples. They show no low frequency localization effects. 



we use l3C 



K(u,n-,T) 



2wcoth(ft/2T) - (w + fi) 



x coth((w + ft)/2T) + (w - ft) coth((w - ft)/2T) 
for the normal state and 
K(u, ft;T = 0) = — ((w - Q.)6(u + 2A (#) - ft) 



(3) 



(4) 



for the superconducting state where (■ ■ • )# denotes the 
angular average over i? and B(x) is the complete elliptic 
integral of the second order. Here A ($) = A cos(2i9) 
reflecting the d-wave symmetry of the superconducting 
order parameter. For the superconducting state at finite 
temperature we adjust the size of the maximum gap, Ao 
according to a BCS temperature variation. After small 
adjustments using full Eliashberg formalism and a least 
squares fit, we obtain the electron-boson spectral density, 
7 2 x(w). Details about uniqueness and quality of the fit 
are described in the literature [sol ]. 

We used the self-energy data of Bi-2212 from a previ- 
ous study jl^ on an optimally doped and two overdoped 
samples. We note that the optimally doped sample is 
yttrium-doped and somewhat different from conventional 
optimally doped Bi-2212 having a more ordered structure 
and a higher T c . Underdoped systems show pseudogap 
behavior not yet incorporated in the inversion formal- 
ism and are not treated here. In Fig. ^ we show our fit 
(dashed lines) to data (solid lines) as an example for our 
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FIG. 2: The bosonic spectra, I 2 \{u), of Bi2Sr2CaCu208+s. 
(a - c) The doping and temperature dependent bosonic func- 
tion for three doping levels (T c = 96 K (OPT), 82 K (OD), 
and 60 K (OD)). At room temperature all samples exhibit a 
broad continuum background shown in panel e. On lowering 
the temperature the broad background peak evolves into a 
broad peak with a deep valley above it. The spectral weight 
gained in the peak is roughly balanced by the spectral weight 
loss in the valley (see Fig. ^p) . With increasing doping the 
intensity of the peak and valley are significantly weakened 
(sec Fig- ISJi) • Panel d shows the temperature dependence of 
the frequency of the peak maximum in OPT96A marked by a 
star. Panel f shows the local magnetic susceptibility of under- 
doped Y-123 from neutron scattering 0, While this 
is a different copper oxide from the ones given in panels a, b 
and c we note the qualitative similar temperature evolution 
of panel a. 



optimally doped sample OPT96A at two temperatures: 
heavy lines are for T=300 K (normal) and light lines for 
27 K (superconducting). The fits are very good ex- 
cept for two small spectral regions near zero energy for 
both states and near the overshoot region for the super- 
conducting state. Above 125 meV in the superconduct- 
ing state the theoretical curve is flat while the data show 
a slight depression. The discrepancy near zero energy 
may come from impurity localization which is observed 
in most cuprate systems. We note that the optimally 
doped sample shows the strongest localization but this 
does not affect to the main sharp peak in I 2 \(^)- The 
two overdoped samples show negligible localization (see 
the insets in Fig. 1). 



3 



In Fig. b and c we plot the extracted bosonic spec- 
tra, 7 2 x(w), of Bi-2212. We observe dramatic tempera- 
ture and doping dependencies. The room temperature 
spectra show a broad continuum which exhibits some 
doping dependence as shown in panel e. At lower doping 
levels the broad peak in the spectrum becomes somewhat 
stronger and shifts to lower frequencies, which is consis- 
tent with a previous neutron study [2j on underdoped 
YBa 2 Cu 3 06.6- The temperature and doping dependence, 
as we move from panel a to panel c, is more striking. On 
the lowering of the temperature nearly all the spectral 
weight below 130 meV in the optimally doped material 
is moved to a resonance peak while in the overdoped re- 
gion less and less of the continuum spectral weight is 
transferred to this peak. The peak in the room tempera- 
ture spectrum, at 97 meV in panel a, shifts continuously 
to lower energies as its intensity grows. The panel d 
shows that the peak position is approximately linear in 
temperature above T = 100 K and saturates to 60 meV 
below this temperature. What is notable about this pro- 
cess is that this evolution takes place through a transfer 
of spectral weight from frequencies above the peak leav- 
ing a valley between 100 to 150 meV. In addition, there 
is a well developed 30 meV gap on the low frequency 
side of the peak. Such a gap in the spin susceptibil- 
ity is consistent with the observations of spin-polarized 
neutron scattering of optimally doped Y-123 with T c =91 
K [3^. As shown in Fig. 03 and [2t for higher doping 
levels the gap in the spin susceptibility becomes small 
and the resonance mode can still be discerned although 
it weakens significantly 26] . In our most overdoped sam- 
ple (OD60G) only a small fraction of spectral weight in 
the bosonic spectrum contributes to the temperature re- 
distribution. 

We also show the local (q integrated) magnetic sus- 
ceptibility, Xodd( £ iO' fr° m neutron scattering study of un- 
derdoped Y-123 |5j,|3] in Fig. EF to compare its tempera- 
ture dependent spectral weight redistribution with that 
of the bosonic excitation, I 2 x(u), of our optimally doped 
Bi-2212 (OPT96A). Although they are different copper 
oxide systems and at different doping levels they qualita- 
tively contain common temperature dependent features 
which we described previously. A similar spectral redis- 
tribution of the local magnetic susceptibility has been ob- 
served by Dai et al. ^J. However, it was not widely rec- 
ognized or confirmed by other experiments. Here we find 
that the spectral redistribution occurs mainly in the low 
frequency region below 200 meV and is strongly doping 
dependent on the overdoped side of the phase diagram. 

Alternatively, we can look at the spectrum at low tem- 
perature in panel a of Fig.EJas the development of a large 
130 meV gap in the spin fluctuation spectrum with a 
mode in the middle of the gap at 60 meV. Such a transfor- 
mation has been considered by Abanov and Chubukov Q 
but with the gap arising at 2Ao with Ao the supercon- 
ducting gap. Here we find a much larger gap that starts 
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FIG. 3: T c dependent properties of the peak and the valley 
(see in the text), a) The central peak frequency is propor- 
tional to T c , i.e., Qp e ak = 8.0fesT c , and the dashed line is a 
linearly fitted line, b) The peak and valley are closely con- 
nected, have a very similar T c dependence and vanish at the 
same T c w 50 K. The dashed lines are linearly fitted lines, c) 
The coupling constant, X pea k, vanishes at the same T c as the 
peak and valley. The coupling constant, Xbg, of the continu- 
ous background shows a weaker T c dependence as compared 
to that of the peak. The dash-dotted line is obtained by using 
the linearly fitted center frequency (Fig. [3^) and area under 
the peak (Fig.^b), i.e., X peak ~ 2A peak /Q peak . 



to form already at T = 200 K, well above T c . 

To study quantitatively the changes in the bosonic 
spectra we subtract the room temperature spectrum from 
the lower temperature spectra. We clearly observe a peak 
and a valley in the difference spectra. In Fig. [31 we dis- 
play the T c dependence of the peak and valley in the 
bosonic spectra. Fig. shows the center frequency of 
the peak, ^l pea k versus T c at our lowest measured tem- 
peratures. We find flpeak ~ 8.0ksT c , which is similar to 
the temperature dependence of the magnetic resonance 
mode (f2 res w 5.3/cbT c ). In Fig. 03 we plot the the area 
under the peak and of the valley above it as a function of 
the critical temperature for our lowest temperature (for 
the optimally doped we also include the 30 meV gap as 
part of the valley). The two areas show the same trend; 
both decrease rapidly as T c decreases in the overdoped 
region and their extrapolation show that they vanish at 
the same T c near 50 K. The total spectral weight up to 
400 meV is conserved to within 15 % which indicates 
a small frequency dependence of the coupling constant 
I 2 in I 2 x(ty- We also calculate the T c dependent cou- 
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pling constant, A(T C ) = 2 f Q c I 2 x(£i) /SldQ, using only 
the peak (X pea k) or the continuum (Xbg) for I 2 x(Q)- 
Here uj c is a cutoff frequency, taken as 400 meV. These 
are shown in Fig. The peak coupling constant X pea k 
decreases rapidly extrapolating to zero at T c =50 K. The 
coupling constant of the background Xbg shows a weaker 
dependence on T c . 

Our novel analysis of the optical data has provided new 
and detailed information on the redistribution with tem- 
perature and doping levels of the spectral weight in the 
inelastic scattering spectral function. This analysis of op- 
tical data, based on an Eliashberg inversion, is the only 
existing method to date which can be used to study the 
evolution of the spin susceptibility as a function of dop- 
ing and temperature across the phase diagram of high-T c 
oxides. Furthermore this analysis can also be applied 
to other existing optical spectra to get new information 
on their excitation spectrum. Our central new finding is 
that the optical resonance peak which becomes promi- 
nent at low temperature in our optimally doped sam- 
ple, forms mainly through the transfer of spectral weight 
from the energy region immediately above it. This con- 
stitute strong evidence that the peak and background in 
the bosonic spectral function have the same microscopic 
origin. Identifying the optical peak with the resonance 
in the local magnetic susceptibility measured by neutron 
scattering, we conclude that the background which domi- 
nates the inelastic interaction at and above T c is magnetic 
in origin. 

Finally we note that we have also identified a new tem- 
perature scale, around 200 K in the optimally doped Bi- 
2212, where the magnetic fluctuation spectrum begins to 
soften and develop the peak in the local magnetic sus- 
ceptibility. To our knowledge current theories do not 
account for this phenomenon. 
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